Quantifying and Cataloguing Unknown Sequences within Human Microbiomes.

Advances in genome sequencing technologies and lower costs have enabled the exploration of a multitude of known and novel environments and microbiomes. This has led to an exponential growth in the raw sequence data that are deposited in online repositories. Metagenomic and metatranscriptomic data sets are typically analysed with regard to a specific biological question. However, it is widely acknowledged that these data sets are comprised of a proportion of sequences that bear no similarity to any currently known biological sequence, and this so-called "dark matter" is often excluded from downstream analyses. In this study, a systematic framework was developed to assemble, identify, and measure the proportion of unknown sequences present in distinct human microbiomes. This framework was applied to 40 distinct studies, comprising 963 samples, and covering 10 different human microbiomes including fecal, oral, lung, skin, and circulatory system microbiomes. We found that while the human microbiome is one of the most extensively studied, on average 2% of assembled sequences have not yet been taxonomically defined. However, this proportion varied extensively among different microbiomes and was as high as 25% for skin and oral microbiomes that have more interactions with the environment. A rate of taxonomic characterization of 1.64% of unknown sequences being characterized per month was calculated from these taxonomically unknown sequences discovered in this study. A cross-study comparison led to the identification of similar unknown sequences in different samples and/or microbiomes. Both our computational framework and the novel unknown sequences produced are publicly available for future cross-referencing. Our approach led to the discovery of several novel viral genomes that bear no similarity to sequences in the public databases. Some of these are widespread as they have been found in different microbiomes and studies. Hence, our study illustrates how the systematic characterization of unknown sequences can help the discovery of novel microbes, and we call on the research community to systematically collate and share the unknown sequences from metagenomic studies to increase the rate at which the unknown sequence space can be classified.

NCBI's Virus Discovery Codeathon: Building "FIVE" -The Federated Index of Viral Experiments API Index.

Viruses represent important test cases for data federation due to their genome size and the rapid increase in sequence data in publicly available databases. However, some consequences of previously decentralized (unfederated) data are lack of consensus or comparisons between feature annotations. Unifying or displaying alternative annotations should be a priority both for communities with robust entry representation and for nascent communities with burgeoning data sources. To this end, during this three-day continuation of the Virus Hunting Toolkit codeathon series (VHT-2), a new integrated and federated viral index was elaborated. This Federated Index of Viral Experiments (FIVE) integrates pre-existing and novel functional and taxonomy annotations and virus-host pairings. Variability in the context of viral genomic diversity is often overlooked in virus databases. As a proof-of-concept, FIVE was the first attempt to include viral genome variation for HIV, the most well-studied human pathogen, through viral genome diversity graphs. As per the publication of this manuscript, FIVE is the first implementation of a virus-specific federated index of such scope. FIVE is coded in BigQuery for optimal access of large quantities of data and is publicly accessible. Many projects of database or index federation fail to provide easier alternatives to access or query information. To this end, a Python API query system was developed to enhance the accessibility of FIVE.

Pingu virus: A new picornavirus in penguins from Antarctica.

Picornaviridae family comprises single-stranded, positive-sense RNA viruses distributed into forty-seven genera. Picornaviruses have a broad host range and geographic distribution in all continents. In this study, we applied a high-throughput sequencing approach to examine the presence of picornaviruses in penguins from King George Island, Antarctica. We discovered and characterized a novel picornavirus from cloacal swab samples of gentoo penguins (Pygoscelis papua), which we tentatively named Pingu virus. Also, using RT-PCR we detected this virus in 12.9 per cent of cloacal swabs derived from P. papua, but not in samples from adélie penguins (Pygoscelis adeliae) or chinstrap penguins (Pygoscelis antarcticus). Attempts to isolate the virus in a chicken cell line and in embryonated chicken eggs were unsuccessful. Our results expand the viral diversity, host range, and geographical distribution of the Picornaviridae.

Metaviromics Reveals Unknown Viral Diversity in the Biting Midge Culicoides impunctatus.

Biting midges (Culicoides species) are vectors of arboviruses and were responsible for the emergence and spread of Schmallenberg virus (SBV) in Europe in 2011 and are likely to be involved in the emergence of other arboviruses in Europe. Improved surveillance and better understanding of risks require a better understanding of the circulating viral diversity in these biting insects. In this study, we expand the sequence space of RNA viruses by identifying a number of novel RNA viruses from Culicoides impunctatus (biting midge) using a meta-transcriptomic approach. A novel metaviromic pipeline called MetaViC was developed specifically to identify novel virus sequence signatures from high throughput sequencing (HTS) datasets in the absence of a known host genome. MetaViC is a protein centric pipeline that looks for specific protein signatures in the reads and contigs generated as part of the pipeline. Several novel viruses, including an alphanodavirus with both segments, a novel relative of the Hubei sobemo-like virus 49, two rhabdo-like viruses and a chuvirus, were identified in the Scottish midge samples. The newly identified viruses were found to be phylogenetically distinct to those previous known. These findings expand our current knowledge of viral diversity in arthropods and especially in these understudied disease vectors.

Krykféie dicistrovirus: A novel dicistrovirus in velvety free-tailed bats from Brazil.

The Dicistroviridae family comprises positive single-stranded RNA viruses that are classified into Picornavirales order. These viruses are identified in arthropod hosts, including some having devastating economic consequences. Here, we described and characterized a novel nearly complete dicistrovirus genome identified in liver samples of velvety free-tailed bats (Molossus molossus) collected in June 2010 in Araçatuba city, São Paulo State, Brazil. This novel virus presents a genome of 9262 nucleotides in length and a typical dicistrovirus genome organization. Based on our phylogenetic analysis and ICTV criteria, we propose this virus as a novel species into the Triatovirus genus. Attempts of viral propagation in Vero E6 and C6/36 cell lines were unsuccessful. The novel dicistrovirus was detected only in one out of nine liver bat samples, representing for the first time an internal organ detection from a representative of this virus family.

Analysis of Paramyxovirus Transcription and Replication by High-Throughput Sequencing.

We have developed a high-throughput sequencing (HTS) workflow for investigating paramyxovirus transcription and replication. We show that sequencing of oligo(dT)-selected polyadenylated mRNAs, without considering the orientation of the RNAs from which they had been generated, cannot accurately be used to analyze the abundance of viral mRNAs because genomic RNA copurifies with the viral mRNAs. The best method is directional sequencing of infected cell RNA that has physically been depleted of ribosomal and mitochondrial RNA followed by bioinformatic steps to differentiate data originating from genomes from viral mRNAs and antigenomes. This approach has the advantage that the abundance of viral mRNA (and antigenomes) and genomes can be analyzed and quantified from the same data. We investigated the kinetics of viral transcription and replication during infection of A549 cells with parainfluenza virus type 2 (PIV2), PIV3, PIV5, or mumps virus and determined the abundances of individual viral mRNAs and readthrough mRNAs. We found that the mRNA abundance gradients differed significantly between all four viruses but that for each virus the pattern remained relatively stable throughout infection. We suggest that rapid degradation of non-poly(A) mRNAs may be primarily responsible for the shape of the mRNA abundance gradient in parainfluenza virus 3, whereas a combination of this factor and disengagement of RNA polymerase at intergenic sequences, particularly those at the NP:P and P:M gene boundaries, may be responsible in the other viruses.IMPORTANCE High-throughput sequencing (HTS) of virus-infected cells can be used to study in great detail the patterns of virus transcription and replication. For paramyxoviruses, and by analogy for all other negative-strand RNA viruses, we show that directional sequencing must be used to distinguish between genomic RNA and mRNA/antigenomic RNA because significant amounts of genomic RNA copurify with poly(A)-selected mRNA. We found that the best method is directional sequencing of total cell RNA, after the physical removal of rRNA (and mitochondrial RNA), because quantitative information on the abundance of both genomic RNA and mRNA/antigenomes can be simultaneously derived. Using this approach, we revealed new details of the kinetics of virus transcription and replication for parainfluenza virus (PIV) type 2, PIV3, PIV5, and mumps virus, as well as on the relative abundance of the individual viral mRNAs.

The Third Annual Meeting of the European Virus Bioinformatics Center.

The Third Annual Meeting of the European Virus Bioinformatics Center (EVBC) took place in Glasgow, United Kingdom, 28-29 March 2019. Virus bioinformatics has become central to virology research, and advances in bioinformatics have led to improved approaches to investigate viral infections and outbreaks, being successfully used to detect, control, and treat infections of humans and animals. This active field of research has attracted approximately 110 experts in virology and bioinformatics/computational biology from Europe and other parts of the world to attend the two-day meeting in Glasgow to increase scientific exchange between laboratory- and computer-based researchers. The meeting was held at the McIntyre Building of the University of Glasgow; a perfect location, as it was originally built to be a place for "rubbing your brains with those of other people", as Rector Stanley Baldwin described it. The goal of the meeting was to provide a meaningful and interactive scientific environment to promote discussion and collaboration and to inspire and suggest new research directions and questions. The meeting featured eight invited and twelve contributed talks, on the four main topics: (1) systems virology, (2) virus-host interactions and the virome, (3) virus classification and evolution and (4) epidemiology, surveillance and evolution. Further, the meeting featured 34 oral poster presentations, all of which focused on specific areas of virus bioinformatics. This report summarizes the main research findings and highlights presented at the meeting.

Discovery of novel astrovirus and calicivirus identified in ruddy turnstones in Brazil.

Birds are the natural reservoir of viruses with zoonotic potential, as well as contributing to the evolution, emergence, and dissemination of novel viruses. In this study, we applied a high-throughput screening approach to identify the diversity of viruses in 118 samples of birds captured between October 2006 to October 2010 in the North and Northeast regions of Brazil. We found nearly complete genomes of novel species of astrovirus and calicivirus in cloacal swabs of ruddy turnstones (Arenaria interpres) collected in Coroa do Avião islet, Pernambuco State. These viruses are positive-sense single-stranded RNA with a genome of ~7 to 8 kb, and were designated as Ruddy turnstone astrovirus (RtAstV) and Ruddy turnstone calicivirus (RTCV), respectively. Phylogenetic analysis showed that RtAstV and RTCV grouped in a monophyletic clade with viruses identified from poultry samples (i.e., chicken, goose, and turkey), including viruses associated with acute nephritis in chickens. Attempts of viral propagation in monkey and chicken cell lines for both viruses were unsuccessful. Also, we found genomes related with viral families that infect invertebrates and plants, suggesting that they might be ingested in the birds' diet. In sum, these findings shed new light on the diversity of viruses in migratory birds with the notable characterization of a novel astrovirus and calicivirus.

A Novel Hepacivirus in Wild Rodents from South America.

The Hepacivirus genus comprises single-stranded positive-sense RNA viruses within the family Flaviviridae. Several hepaciviruses have been identified in different mammals, including multiple rodent species in Africa, Asia, Europe, and North America. To date, no rodent hepacivirus has been identified in the South American continent. Here, we describe an unknown hepacivirus discovered during a metagenomic screen in Akodon montensis, Calomys tener, Oligoryzomys nigripes, Necromys lasiurus, and Mus musculus from São Paulo State, Brazil. Molecular detection of this novel hepacivirus by RT-PCR showed a frequency of 11.11% (2/18) in Oligoryzomys nigripes. This is the first identification of hepavivirus in sigmondonine rodents and in rodents from South America. In sum, our results expand the host range, viral diversity, and geographical distribution of the Hepacivirus genus.

ICTV Virus Taxonomy Profile: Parvoviridae.

Members of the family Parvoviridae are small, resilient, non-enveloped viruses with linear, single-stranded DNA genomes of 4-6 kb. Viruses in two subfamilies, the Parvovirinae and Densovirinae, are distinguished primarily by their respective ability to infect vertebrates (including humans) versus invertebrates. Being genetically limited, most parvoviruses require actively dividing host cells and are host and/or tissue specific. Some cause diseases, which range from subclinical to lethal. A few require co-infection with helper viruses from other families. This is a summary of the International Committee on Taxonomy of Viruses (ICTV) Report on the Parvoviridae, which is available at www.ictv.global/report/parvoviridae.

Viral diversity of Rhipicephalus microplus parasitizing cattle in southern Brazil.

Ticks are ectoparasites spread worldwide and are well known as vectors of many viruses of great importance to human and animal health. However, the viral diversity in ticks is still poorly understood, particularly in South America. Here we characterized the viral diversity present in Rhipicephalus microplus parasitizing cattle in the southern region of Brazil using metagenomics. Our study revealed the presence of viruses that had not been previously described in the region, including lihan tick virus (Phenuiviridae family) and wuhan tick virus 2 (Chuviridae family), as well as expands the biogeography of jingmen tick virus (Flaviviridae family) in Brazil. Also, we described three novel tymoviruses (Tymovirales order), named guarapuava tymovirus-like 1 to 3. We described the genomic and phylogenetic characterization of these viruses. Our study sheds light on the viral diversity of Rhipicephalus microplus in South America, and also expands the biogeography of tick viruses that were previously described only in Asia.

Development of reverse genetics systems and investigation of host response antagonism and reassortment potential for Cache Valley and Kairi viruses, two emerging orthobunyaviruses of the Americas.

Orthobunyaviruses such as Cache Valley virus (CVV) and Kairi virus (KRIV) are important animal pathogens. Periodic outbreaks of CVV have resulted in the significant loss of lambs on North American farms, whilst KRIV has mainly been detected in South and Central America with little overlap in geographical range. Vaccines or treatments for these viruses are unavailable. One approach to develop novel vaccine candidates is based on the use of reverse genetics to produce attenuated viruses that elicit immune responses but cannot revert to full virulence. The full genomes of both viruses were sequenced to obtain up to date genome sequence information. Following sequencing, minigenome systems and reverse genetics systems for both CVV and KRIV were developed. Both CVV and KRIV showed a wide in vitro cell host range, with BHK-21 cells a suitable host cell line for virus propagation and titration. To develop attenuated viruses, the open reading frames of the NSs proteins were disrupted. The recombinant viruses with no NSs protein expression induced the production of type I interferon (IFN), indicating that for both viruses NSs functions as an IFN antagonist and that such attenuated viruses could form the basis for attenuated viral vaccines. To assess the potential for reassortment between CVV and KRIV, which could be relevant during vaccination campaigns in areas of overlap, we attempted to produce M segment reassortants by reverse genetics. We were unable to obtain such viruses, suggesting that it is an unlikely event.

A novel polyomavirus in sigmodontine rodents from São Paulo State, Brazil.

The nearly complete genome sequence of a novel polyomavirus from blood samples of Akodon montensis and Calomys tener collected in Brazil was determined by high-throughput sequencing. This virus showed a typical polyomaviruses genome organization, and it was classified as a member of the genus Betapolyomavirus. Our results expand the host range and viral diversity of the family Polyomaviridae.

Novel orthohepeviruses in wild rodents from São Paulo State, Brazil.

The Hepeviridae comprise single-stranded positive-sense RNA viruses classified into two genera, Orthohepevirus and Piscihepevirus. Orthohepeviruses have a wide host range that includes rodents, but previous studies had been restricted to rodents of the Muridae family. In this study, we applied a high-throughput sequencing approach to examine the presence of orthohepeviruses in rodents from São Paulo State, Brazil. We also used RT-PCR to determine the frequency of orthohepeviruses in our sampled population. We identified novel orthohepeviruses in blood samples derived from Necromys lasiurus (1.19%) and Calomys tener (3.66%). Therefore, our results expand the host range and viral diversity of the Hepeviridae family.

Novel Parvoviruses from Wild and Domestic Animals in Brazil Provide New Insights into Parvovirus Distribution and Diversity.

Parvoviruses (family Parvoviridae) are small, single-stranded DNA viruses. Many parvoviral pathogens of medical, veterinary and ecological importance have been identified. In this study, we used high-throughput sequencing (HTS) to investigate the diversity of parvoviruses infecting wild and domestic animals in Brazil. We identified 21 parvovirus sequences (including twelve nearly complete genomes and nine partial genomes) in samples derived from rodents, bats, opossums, birds and cattle in Pernambuco, São Paulo, Paraná and Rio Grande do Sul states. These sequences were investigated using phylogenetic and distance-based approaches and were thereby classified into eight parvovirus species (six of which have not been described previously), representing six distinct genera in the subfamily Parvovirinae. Our findings extend the known biogeographic range of previously characterized parvovirus species and the known host range of three parvovirus genera (Dependovirus, Aveparvovirus and Tetraparvovirus). Moreover, our investigation provides a window into the ecological dynamics of parvovirus infections in vertebrates, revealing that many parvovirus genera contain well-defined sub-lineages that circulate widely throughout the world within particular taxonomic groups of hosts.

ViCTree: an automated framework for taxonomic classification from protein sequences.

Motivation: The increasing rate of submission of genetic sequences into public databases is providing a growing resource for classifying the organisms that these sequences represent. To aid viral classification, we have developed ViCTree, which automatically integrates the relevant sets of sequences in NCBI GenBank and transforms them into an interactive maximum likelihood phylogenetic tree that can be updated automatically. ViCTree incorporates ViCTreeView, which is a JavaScript-based visualization tool that enables the tree to be explored interactively in the context of pairwise distance data. Results: To demonstrate utility, ViCTree was applied to subfamily Densovirinae of family Parvoviridae. This led to the identification of six new species of insect virus. Availability and implementation: ViCTree is open-source and can be run on any Linux- or Unix-based computer or cluster. A tutorial, the documentation and the source code are available under a GPL3 license, and can be accessed at http://bioinformatics.cvr.ac.uk/victree_web/. Supplementary information: Supplementary data are available at Bioinformatics online.

Discovery of novel anelloviruses in small mammals expands the host range and diversity of the Anelloviridae.

The Anelloviridae comprises single-stranded DNA viruses currently grouped in sixty-eight species classified in twelve genera. They have been found in many vertebrate hosts including primates. In this study, we describe the application of the high-throughput sequencing to examine the frequency and diversity of anelloviruses in rodents, bats and opossums captured in São Paulo State, Brazil. We report a total of twenty-six anelloviruses with sixteen nearly complete genomes and ten partial genomes, which include eleven potential novel species identified in rodents (Cricetidae), bats (Molossidae and Phyllostomidae), and opossums (Didelphidae). We also propose the inclusion of two potential new genera within the Anelloviridae family, provisionally named Omegatorquevirus and Sigmatorquevirus, including six and three novel species of anelloviruses, respectively. In summary, this study expands the diversity and the host range of the known anelloviruses.

Chapparvoviruses occur in at least three vertebrate classes and have a broad biogeographic distribution.

Chapparvoviruses are a highly divergent group of parvoviruses (family Parvoviridae) that have recently been identified via metagenomic sampling of animal faeces. Here, we report the sequences of six novel chapparvoviruses identified through both metagenomic sampling of bat tissues and in silico screening of published vertebrate genome assemblies. The novel chapparvoviruses share several distinctive genomic features and group together as a robustly supported monophyletic clade in phylogenetic trees. Our data indicate that chapparvoviruses have a broad host range in vertebrates and a global distribution.

The use of next generation sequencing in the diagnosis and typing of respiratory infections.

BACKGROUND: Molecular assays are the gold standard methods used to diagnose viral respiratory pathogens. Pitfalls associated with this technique include limits to the number of targeted pathogens, the requirement for continuous monitoring to ensure sensitivity/specificity is maintained and the need to evolve to include emerging pathogens. Introducing target independent next generation sequencing (NGS) could resolve these issues and revolutionise respiratory viral diagnostics. OBJECTIVES: To compare the sensitivity and specificity of target independent NGS against the current standard diagnostic test. STUDY DESIGN: Diagnostic RT-PCR of clinical samples was carried out in parallel with target independent NGS. NGS sequences were analyzed to determine the proportion with viral origin and consensus sequences were used to establish viral genotypes and serotypes where applicable. RESULTS: 89 nasopharyngeal swabs were tested. A viral pathogen was detected in 43% of samples by NGS and 54% by RT-PCR. All NGS viral detections were confirmed by RT-PCR. CONCLUSIONS: Target independent NGS can detect viral pathogens in clinical samples. Where viruses were detected by RT-PCR alone the Ct value was higher than those detected by both assays, suggesting an NGS detection cut-off - Ct=32. The sensitivity and specificity of NGS compared with RT-PCR was 78% and 80% respectively. This is lower than current diagnostic assays but NGS provided full genome sequences in some cases, allowing determination of viral subtype and serotype. Sequencing technology is improving rapidly and it is likely that within a short period of time sequencing depth will increase in-turn improving test sensitivity.